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INTRODUCTION 



The importance of the double-well potentials in quantum mechanics, condensed matter physics, statistical physics 
or field theory can hardly be overestimated [1, 2]. The Schrodinger equation has no exact, analytical solution with 
double-well potentials. It can be solved by using approximate methods or numerical methods. On the other hand, 
in quantum mechanics there exist potentials for which it is possible to find a number of eigenvalues and associated 
eigenfunctions exactly and in a closed form. These systems are said to be quasi-exactly solvable (QES) and this 
property is ultimately connected with the existence of a hidden dynamical algebra [3] . 

In this article we present a quasi exact solution of the Schrodinger equation with one- and two dimensional double- 
well potentials. The Schrodinger equation with one dimensional double- well potential have been studied by a number 
of authors and it can be solved by using various methods [4] . Approximate QES and approximate analytical solutions 
of the double- well potential have been derived in [5] . 

There is not much attention for the solution of the two-dimensional double- well potential [6]. In this paper we 
develop a simple method to solve the Schrodinger equation which admits separation of variables with two-dimensional 
double-well potentials. We show the separated equations are still QES. The Schrodinger equation does not admit 
separation of variables in two dimensions, a topic of a future work. 

The paper is organized as follows. Section 2 is devoted to a derivation of a QES one-dimensional double-well 
potential in the framework of the SI2 {R) Lie algebra. Section 3 contains a solution of the two dimensional double- well 
potential which admits separation of variables. Finally, the results are discussed in section 4. 

ONE DIMENSIONAL DOUBLE- WELL POTENTIAL 

The strategy we follow is described in [7] . One way to construct a quasi-exactly solvable differential equation is to 
express the Hamiltonian in terms of the generators of the Lie algebra. Let us consider the SI2 (R) algebra realized as 
follows: 

2d d j d 

The generators obey the commutation relation 

[J+, J_] = -2Jo, [Jo, ^±]=±J±, (2) 
and it may be verified that the eigensolutions of the generators are given explicitly by 

J_| m >= m I m — 1 >, Jo | to ( to — — ) | m>, J+ | m>= (to — j)\m + 1), (3) 
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where the basis function is 



I TO >= Rj+i{x). 



(4) 
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If j is a non-negative integer number, the algebra possesses a (j + l)-diniensional subspace 

Rj+i = {x°,x'',--- ,x^)- (5) 

The linear and bilinear combinations of the generators of the algebra form a second order differential equation 

possessing polynomial solutions. One possible way to get a QES Schrodinger operator is to transform the QES algebraic 
operators acting on the finite dimensional subspace of polynomials Kj-^-l{x) into the Schrodinger-type equation. It 
can be done by a gauge transformation. 



The model Potential I 



Consider the double-well potential of the form 

where A, B and C arc constants. If 2B^ = 9AC, the potential is a symmetric potential with two minima; otherwise 
the potential is an asymmetric potential with two minima. We are interested in the asymmetric double-well potential. 
The Schrodinger equation of this potential can be written as 

.^ + V{x)^^{x)=E^{x). (7) 

We can write the wavefunction 'tjj{x) in the form 

Hx) = exp (^-30 + 1)_ - + j Rix). (8) 

Then the Schrodinger equation takes the form 

r-^ + -3(i + l)-^ i=x + \/Cx^ — i-^ - 

dx^ V B 3\/C ) dx 

with the constraint 



^ VCx ] ] R{x) = (9) 



A = 



{2B^ - 27(j + 1)CV2) 



9BC 
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E = e-il±^-('-^l±^) . (10) 



3VC V 2B 

The potential (6) is asymmetric for all values of j when B > and C > 0, and the location of minimum points 
depends on the values of j. In order to solve (9), we introduce the linear and bilinear combination of the algebraic 
operators 

T = -Jl- 3(i + 1)- J_ - ^ Jo + VCJ+ (11) 
for which one can define the spectral problem 

TR(a;) = eR{x), (12) 

where e is a spectral parameter. The algebraic structure (11) is quasi-exactly solvable and the insertion of (1) into 
(11) leads to the differential equation (9). The basis function (5) can be written as 

oo 

R{x) = amPm{e)x"', (13) 

m=0 
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where Pmi^) is polynomial in e and it can be given by the recurrence relation 

C 

-m{m- 1)P™_2(£) - 3rn(j + l)—Pm-iie) + Vc {m - j) Pm+i{e) - 

(^l^{m-l)+e^P^{e) = Q, (14) 

with the initial condition Pq = 1. The spectral parameter e can be determined from the recurrence relation (14). If 
Ej is a root of the polynomial P„i+i{e), the scries (13) truncates at m ^ (j + 1) and Ej belongs to the spectrum of 
the double-well potential. The first four of these polynomials are given by 

Pi{e) = E 

Pz{e) = gSCe^ - 4 + 81C^/2 j e + 36SC2 
Pi{£) = QS^CV- (105^(7+ 1080BC^/2) 

216B^C^E + (^B^ + 216B^C^/^ + 11664C^) (15) 

for j = 0, 1, 2 and 3, respectively. Analytical solutions of the recurrence relation (14) are available only for the first 
few values of j. For j > 2 the solutions become numerical. 

The Model Potential II 

We consider the following potential treated by Burrows [5] 

V{x)=Vo-Ax + —{B + Cxf. (16) 

The potential is symmetric when A = 0, with double minima at a; = — ^ and x = The Schrodinger equation 

with the potential (16) is given by 

^ + V{x)^ i>{x) = E^ix). (17) 

In order to transform the Schrodinger equation to the form of an algebraic equation, we introduce the following 
transformation: 

i;{x) = exp (-^ - i?(x). (18) 



2V2 3V2 



Then the Schrodinger equation takes the form 



-'^-V2xiB + Cx)'-^ + 
dx-^ ax 

(vq-E-Ax + V2Cx - R{x) = 0. 



(19) 



The linear and bilinear combinations of the operators of the Lie algebra 

T =-Jl~ V2BJ0 - V2CJ+; TR{x) = eR{x) (20) 
is equivalent to the situation in (19) when the parameters of the Schrodinger equation is constrained to 

A = -s/2C{j + l) 

E = E-^1±^. (21) 
^/2 
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Since j ^ 0, the potential is an asymmetric double-well potential. The Spectral parameter e can be calculated from 
the recurrence relation 

-m(m - l)P™-2(£) - (^V2B (m - + Pm{e) - (22) 

V2C{m-j)Pm+i{e) = 0, 
with the initial condition Po{£) = 1- The first four values of the Pm{s:) are given as follows: 

Pi{e) = e 

P2{e) = 2e^-B^ 

P3(e) = V2{e^ -2eB^ -ABC) +4:eC 

Piie) = As" - 20 - 2V2C) + 96£BC + 95^ _ ^^c^ , (23) 
for j = 0, 1, 2, 3, respectively. 



TWO-DIMENSIONAL DOUBLE- WELL POTENTIAL 



In this section the quasi exact solution of the Schrodinger equation with two dimensional double-well potential is 
treated. We develop a systematic procedure for constructing a QES of a double-well potential and we illustrate our 
method in the situation where the Schrodinger equation admits separation of variables. The Schrodinger equation in 
two dimensions is given by 

-^-^^V{x,y)\^^{x,y)=Ei,{x,y), (24) 

and can be solved exactly for a few potentials. The ground state wave function of the Schrodinger equation can be 
written in the form 

^o(a;, V) = exp(-W^(ar) - F{y) - f{xy)), (25) 

where W{x) and F{y) are functions of x and y, respectively, while f{xy) is function of the product of x and y. The 
potential and ground state energy of the system with the given ground state wave function takes the form 

dx^ dy"^ dx? dy^ 

When f{xy) is zero, then the equation is separable and the general solution of the system can be obtained by solving 
two one dimensional differential equation. The wave function of the excited states can be written as 

il){x,y) = (t){x,y)'il)oix,y), (27) 

where the function ^{x, y) is polynomial in x and y. Substituting (27) into Schrodinger equation (24) we obtain the 
following differential equation: 



Voix,y)-Eo = —^ + -T + w^ + wi-{-^ + ^] + ■ (26) 



(dW^ 








\ dx 


dx ) 


ax \ ay 


dyj 



dx^ dy^ \ dx ' dx J dx \dy dy J dy~^ ^ 
{V -Vo + Eq-E)<I> = Q. 

In the exactly solvable case, the potentials V and Vq are same. In the quasi exactly solvable case, the functional form 

of the potential remain same for all states, but a parameter changes from state to state. Let us illustrate our method 
on the separable systems. In the case / = the potentials can be written as 

V -Vo=9{x)+h{y), (29) 

then substituting = R(x)Q{y) we obtain 

d'^R dW dR , , ^ ^ 

d^Q dFdQ , 
In the following we illustrate our method on two examples. 
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The model potential I 

Consider the following model potential with the parameter a: 

V{x, y; a) = A^x^ + ABx^ + C^j/ + 2Aax. (31) 

The Schrodinger equation with this potential is separable and can be solved by reduction to two one-dimensional 
systems. The parameter a is depends on the state of the wave function. Let us call it the "state parameter". The 
wave function of the ground state with the energy Eq can be written in the form 

Wo(x,y) = exp I Bx 

Eo = C-B"^. (32) 
The wave functions of the excited states with the energy E we write as 

i}{x,y) = R{x)Q{y)Mx,y)- (33) 

Substituting (30) in (24) with the potential (31), then after a straightforward calculation we obtain two one- 
dimensional differential equations: 

-'^+2{C + Ax'')'^ + {E-Eo + c^+2Ax{l-a))R = 0, (35) 
where ci is a constant. The first equation is exactly solvable and its solution is given by 

Q{y)=N,,F, Q-|L,|fiy2^ + N,,F, , (36) 

where Ni and N2 are normalization constants and iFi{---) is a confluent hypcrgcomctric fimction. The second 
equation is a quasi-exactly solvable differential equation. The following combinations of the operators of the Lie 
algebra 

T = -j'i + 2BJ_ + 2AJ+ TR={Eq-E-cx)R (37) 

is equivalent with (35) if a = j + 1. Following the procedure given in section 2, the spectrum of the system can be 
obtained by using the following recurrence relation: 

-m(m - 1)P^_2{E) + 2BmPm-i{E) + (38) 
2A{m - j)Pm+i{E) -{Eo-E- c^)Pm{E) = 0. 

Some of the polynomial is given by 

Pi{E) = (ci-Eo + E) 

P2{E) = {ci - Eo + Ef + AAB 

PsiE) = (ci -Eo + Ef + 16 AB (ci - Eq + E) - 16 A^ 

Pi{E) = [ci-Ea + Ef + AQAB{cx-Eo + Ef - (39) 
96^2 (ci -Eo + E) + lUA^B'^, 

for j = 0, 1, 2, 3, respectively. The roots of the polynomials Pj{E) produce the spectrum of the system. 



The model potential II 



The last example is the following double-well potential: 

V{x, y; a; /3) = A^x^ + Aly^ + 2AiBix^ + 2^2 Say^ - 2Aiax - 2^2/3?/. 



(40) 



As we mentioned before, a and /3 are state parameters. The ground state wave function is given by 

, . ^ [ Mx^ A^y^ „ \ 
i^o{x,y) = expl Bix-B2yj; 

Eo = Bf + Bl (41) 

Since the Schrodinger equation is separable with the given potential, the wave function can be written as 

ij{x,y)=R{x)Q{y)i;o. (42) 

Insertion of (42) into (24) with the potential (40) leads to the following two differential equations: 

-y^ + 2{Bi+Aix'')^ + {E-Eo + c,+2Ai{l-a)x) = (43) 



dx'^ ^ ' dx 

2 (52 + A2t/')^-(ci- 2^2(1-/3) 2/) = 0. (44) 



.— + 2{B, + A,y) — 



These equations are QES. 



CONCLUSION 



The QES of the double-well potential in one and two dimensions have been discussed. The solution of Schrodinger 
equations which does not admit separation of variables are left for treatment in future work. 
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